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Executive Summary 


The inference phase of rocket engine health monitoring involves analysis of results from the 
acquisition phase, comparison of analysis results to establish health measures, and assessment of 
health indications. A particular model based procedure developed for the inference phase of 
engine health monitoring and referred to as Generalized Data Reduction (GDR) was the subject 
of this effort. The GDR method can be considered a strategy for solving the inverse performance 
analysis problem often referred to as data reduction. The primary objective of this research 
effort was to investigate methods of enhancing the GDR strategy in order to make it a flight 
capable real-time diagnostic platform for the engine system of choice. 

A brief conceptual description of the basic GDR strategy punctuated by significant 
mathematical results and a description of major components is presented. The MC-1 engine 
system was specified as the test platform for evaluating GDR capabilities and assessing 
candidate enhancement procedures. Data from the R2 and R3 series of MC-1 engine tests 
conducted at Rocketdyne’s Santa Susanna Field Laboratory in California were utilized to seed 
the data reduction process. A ROCETS performance model of the MC-1 engine was used to 
generate influence matrices used by GDR and a modified ROCETS data reduction model 
provided conventional data reduction results used to assess the accuracy potential of GDR 
predictions. 

The basic GDR procedure was shown to provide consistent approximation of hardware 
function that agreed well with the recognized standard for MC-1 engine data reduction. 
Modifications to the basic procedure were explored to increase the effective solution range as 
well as improve computational speed, accuracy, and stability of the solution procedure. An 
efficient procedure based on singular value decomposition (SVD) of the hardware influence 
matrix was implemented and tested as was a partial second order extension of the basic GDR 
procedure. Computational efficiency was improved using the SVD technique. No improvement 
was observed using the partial second order extension. 

The original MC-1 hardware parameter set used for data reduction was functionally 
regrouped and GDR analyses were performed with the revised hardware set. Agreement of 
predictions with accepted standards was excellent. Single source anomaly resolution capability 
was found to be excellent using the revised parameter set. A study was performed to determine 
the effects of single sensor failure on reduction predictions. Degradation of reduction results in 
the presence of sensor failure was characterized. GDR analyses were also performed using a 
highly restricted flight measurement suite. Loss of hardware discrimination capability using 
only flight available measurements was characterized and results were evaluated as reasonable 
for real-time monitoring applications. 

Two automated subset selection procedures were developed in order to identify compatible 
sets of measurements and hardware parameters for data reduction applications. Results of 
preliminary testing are reported and the shortcomings of sequential subset selection procedures 
are identified. 

Summary evaluation of GDR performance for MC-1 data reduction is provided and 
recommendations for future development are presented. 
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1.0 Introduction 


Monitoring the health of rocket engine systems is essentially a two-phase process. The 
acquisition phase involves sensing physical conditions at selected locations, converting physical 
inputs to electrical signals, conditioning the signals as appropriate to establish scale or filter 
interference, and recording results in a form that is easy to interpret. The inference phase 
involves analysis of results from the acquisition phase, comparison of analysis results to 
established health measures, and assessment of health indications. 

A variety of analytical tools may be employed in the inference phase of health monitoring. 
These tools can be separated into three broad categories: statistical, rule based, and model based. 
Statistical methods can provide excellent comparative measures of engine operating health. 

They require well-characterized data from an ensemble of “typical” engines, or “golden” data 
from a specific test assumed to define the operating norm, in order to establish reliable 
comparative measures. Statistical methods are generally suitable for real-time health monitoring 
because they do not deal with the physical complexities of engine operation. The utility of 
statistical methods in rocket engine health monitoring is hindered by practical limits on the 
quantity and quality of available data. This is due to the difficulty and high cost of data 
acquisition, the limited number of available test engines, and the problem of simulating flight 
conditions in ground test facilities. In addition, statistical methods incur a penalty for 
disregarding flow complexity and are therefore limited in their ability to define performance shift 
causality. 

Rule based methods infer the health state of the engine system based on comparison of 
individual measurements or combinations of measurements with defined health norms or rules. 
This does not mean that rule based methods are necessarily simple. Although binary yes-no 
health assessment can sometimes be established by relatively simple rules, the causality 
assignment needed for refined health monitoring often requires an exceptionally complex rule 
base involving complicated logical maps. Structuring the rule system to be clear and 
unambiguous can be difficult, and the expert input required to maintain a large logic network and 
associated rule base can be prohibitive. 

Model based methods incorporate physical relations and empiricisms in the inference phase 
of health monitoring. Such methods are typically more involved because the flow physics of 
rocket engines is generally described by complex and interdependent, nonlinear relations. The 
attending computational complexity presents a significant impediment to the use of model based 
methods in real-time health monitoring. However, the addition of physical detail does provide a 
basis for determining performance shift causality at the component level. Recent surveys of 
model based inference procedures applicable to propulsion systems are available [1,2]. 

A particular model based procedure referred to as Generalized Data Reduction (GDR) [3, 4] 
was the subject of this effort. GDR was initially developed at Marshall Space Flight Center 
(MSFC). Mathematically, GDR can be considered a strategy for solving the inverse performance 
analysis problem often referred to as data reduction. More specifically, the GDR method inverts 
a canonical representation of the engine system performance model to estimate parameters that 
characterize the operation of hardware components such as pumps, turbines, injectors, nozzles, 
orifices, valves, etc., consistent with test data. As a propulsion system analysis procedure, the 
basic form of GDR is similar to gas path analysis techniques [see e.g. 5] commonly used for fault 
detection and isolation in air-breathing systems. GDR has been applied to mainstage data 
reduction of both the Space Shuttle Main Engine (SSME) [3] and the MC-1 engine system [4]. 
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In these applications GDR demonstrated computational speeds consistent with real-time 
monitoring requirements. 

2.0 Objectives 

The primary objective of the subject research effort was to investigate methods of enhancing 
the Generalized Data Reduction strategy in order to make it a flight capable real-time diagnostic 
platform for the engine system of choice. The principal areas targeted for enhancement were as 
follows: 

1 . Basic algorithm refinement to improve computational efficiency and reduction 
cycle speed 

2. Extensions to incorporate additional engine system knowledge, especially 
uncertainty estimates 

3. Subset selection refinement to improve stability and sensitivity of hardware- 
measurement compatibility assessment and selection process 

4. Range extension methods to expand the engine operating range of application 

5. Techniques for responding to sensor degradation /failure 

6. Methods of incorporating creep scale transients within main stage engine 
diagnostics. 

Secondary topics of study that depended on the GDR development path were identified in the 
grant supplement as described below: 

7. Methods for improving the gateway for input of statistical parameters to the 
data reduction process 

8. Methods for improving data reduction stability with expanded measurement 
and/or hardware sets 

9. Methods for improving performance model feedback from data reduction 
predictions 

10. Methods for improving test data/model accuracy assessment 

1 1 . Methods for improving reduction model capability with limited measurement 
sets. 

A brief description of the basic GDR method is presented in the next section of this 
report in order to introduce terminology and define the development baseline. 

3.0 The GDR Strategy 

A detailed development of the basic GDR strategy can be found in two recent papers [3,4]. 
A brief conceptual description punctuated by significant mathematical results and a description 
of major component procedures is provided below. 

The mainstage performance model of a rocket engine system can be thought of as a set of 
functional relations F between the following parameter sets. 
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1 . known/measured system control settings and/or input/output conditions C; e.g., 
valve positions, inlet/outlet temperatures, pressures, etc. 

2. nominally fixed parameters H describing component hardware function; e.g., 
efficiencies, discharge coefficients, head coefficients, line resistances, etc., and 

3. internal engine conditions P at various locations; e.g., temperatures, pressures, 
flows, heat transfer rates, thrust, etc. 


The objective of performance analysis is to predict the internal engine conditions P, given known 
or assumed values for the members of H and C. 

If the performance model relations between P, H and C are assumed to be linear near some 
known base state designated by the subscript “o”, then they can be approximated in the following 
form. 


Linearized Performance Model Relations 

AP = J Ho AH + Jo AC 


where 


Jho - 


ppil 


1 

> 
i 

_3Hj_ 

0 

L AHjJ 


Jco - 


AH = H-H„ 


AC = C-C 0 


dPi ~ APj 

L<3C k _o . ACkJo 

AP = P-P„ 


H = [H 1 H 2 ...H n ] T C = [C 1 C 2 ...C L ] T P=[P,P 2 ...P m l 1 


( 1 ) 

( 2 ) 

(3) 

(4) 


J Ho mxn Jacobian matrix of P with respect to H at base state o 

Jo mxL Jacobian matrix of P with respect to C at base state o 

The Jacobian matrices defined in relations (2) above are commonly referred to as first 
order influence matrices. Elements can be approximated using a standard finite difference 
scheme and results of performance model simulation runs at incremented values of each 
hardware and control component about its defined base state value. 

The problem of determining values of the hardware shifts AH so that equation system (1) 
is satisfied, given measurement indicated values for AP and AC, is the linear data reduction 
problem. Unfortunately determination of the hardware parameter solution AH is complicated by 
several factors. In general, the number of hardware parameters Hj is greater than the number of 
available measurements Pj, i.e. n>m. This prevents straightforward inversion of the hardware 
Jacobian to obtain a unique solution to equation system (1). It also presents the problem of 
selecting the correct hardware operating state from an infinite set of candidate hardware states 
for the given observation set. Added to the discrimination problem is another related to 
computational stability. If the hardware Jacobian is ill conditioned, the data reduction solution 
AH may be unstable. Little confidence can be placed in the solution to an ill conditioned system 
even with low levels of measurement noise. 

The basic GDR strategy attempts to address these problems using a two-stage solution 
process. Terminology pertinent to the solution procedure is given below. 
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H s vector subset of n s hardware parameters selected for adjustment in reduction analysis 
P s vector subset of m s measured internal conditions selected for reduction analysis matching 

Jhso m s xn s submatrix of hardware Jacobian Jho corresponding to selected hardware, H s , and 

measured internal conditions, P s 
K2(A) L2 condition number of the indicated matrix A 

W n s xn s diagonal matrix of weighting factors associated with selected hardware H s 
QR-Pi column pivoted QR factorization of selected mxn matrix A, i.e., APi = QR where Q is an 
orthogonal mxm matrix, R is a mxn upper triangular matrix, and Pi is a nxr permutation 
matrix 

Using the above terminology, the basic GDR strategy can be described as follows. 

Basic GDR Strategy 

Process 1. Subset Selection 


Select 

H s cH 

and 

P S CP 


using 

QR-Pi 

QR-Pi 

factorization of 
factorization of 

Jho to select 

Jho T to select 

H s cH 
Ps cP 

such that 

K2(Jhso) 

^ ^2-threshold 




Process 2. Optimized Reduction 

Select AH S to 

X 

minimize AH S W AH S 

subject to AP S = Jhso AH s + J Cso AC S 

The subset selection process provides a method of identifying appropriate combinations of 
hardware parameters and internal measurement variables that limit the solution error bound. The 
L2 condition number of the hardware Jacobian, K2 (Jhso)> affords a measure of the computational 
limits of solution accuracy separate from model and measurement uncertainty effects. This is 
reflected in a general guideline for underdetermined linear systems which estimates the number 
of significant digits lost in the solution process to be equal log(K2). 

The pivoted QR factorization process [see, e.g. 6 , 7 ] sequentially arranges hardware 
parameter columns and internal measurement rows of the influence matrix Jh 0 in descending 
order of independence up to the rank of the matrix. After pivoting, elimination of parameters 
associated with later columns and measurements associated with later rows tends to reduce the 
resulting submatrix condition number, thereby reducing the computational error bound. The 
penalty associated with column elimination is artificial reduction in the hardware parameter state 
space range. The penalty associated with row elimination is loss of state space discrimination. 
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The accuracy trade-offs for underdetermined systems are clearly indicated in a perturbation 
result presented by Golub, et al [6, page 273], 

The optimized reduction process requires solution of a weighted least squares problem 
with equality constraints. Orthogonalization methods based on QR factorization or singular 
value decomposition (SVD) are effective and dependable solution procedures for problems of 
this type [6 7 89 10 ]. However, the basic GDR strategy utilizes Powell's implementation of the 
Goldfarb and Idnani dual quadratic programming algorithm [II ] as the solution procedure. This 
method, although less efficient, provides the flexibility to include inequality constraints as well 
as linear and covariant terms in the optimization merit measure. 

A conceptual view of the GDR strategy is presented in the figure below. 


Figure 1 GDR Conceptual Strategy 



Influence 


Baseline 

Matrices 
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Definition 



4.0 Study Components 

To define the study basis, a description of three specific components is needed: the engine 
system that served as the test platform for evaluation of GDR, the test data source for the 
specified engine system, and the parent engine performance model that was used to derive 
engine system influences. Definition of these components is provided below. 
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4.1 Engine System 

The MC-1 engine system was specified as the test platform for evaluating GDR 
capabilities and assessing candidate enhancement procedures. The MC-1 engine is a 60,000 lb 
thrust, pump-fed liquid fuel rocket engine that was developed by the Marshall Space Flight 
Center (MSFC) [12]. It bums a mixture of RP-1 hydrocarbon fuel and liquid oxygen. Hot gas 
produced by a gas generator is used to power a turbine that rotates an inline turbopump 
assembly. The engine uses five orifices to control engine thrust and mixture ratio. These are the 
gas generator liquid oxidizer orifice, the gas generator RP-1 orifice, the main liquid oxidizer 
orifice, the main RP-1 orifice, and the gas generator nozzle orifice. The MC-1 engine is intended 
to be reusable with the exception of the ablative nozzle. A conceptual schematic displaying 
major components of the engine system is presented in Appendix A, Figure A1 . 

4.2 Test Data Source 

Three series of MC-1 tests were recently conducted on the Alfa-1 test stand at 
Rocketdyne’s Santa Susanna Field Laboratory in California. The R1 series was focused on 
resolving any issues associated with the new test stand. The R2 and R3 series consisted of 
several tests of Engine 3 and Engine 5, respectively. Engines 1 and 2 had previously been tested 
in the horizontal test facility at NASA Stennis Space Center (SSC), and Engine 3 had previously 
been tested at the Propulsion Test Article Facility at SSC. The R2 series used a 1 5 : 1 area ratio 
nozzle and the R3 series used a 30:1 area ratio nozzle. Both series were conducted at test stand 
altitude. 

Both R2 and R3 series data sets were utilized in this study. The R2 series consisted of five 
tests: the first four - R2-1, R2-2, R2-3a and R2-3b - were 24 seconds in duration, and the final 
test - R2-4 - was a full duration test of 1 59 seconds. The R3 series consisted of two tests: 

R3-la and R3-2b, each 24 seconds in duration. 

The tests had varying objectives [13]. R2-1 established an Engine 3 calibration baseline. 
The effects of LOX inlet pressure variations on pump run characteristics were investigated in 
R2-2. R2-3a was intended to evaluate engine calibration in response to an orifice change; 
however, due to a leak in the oxidizer bleed valve, this was not accomplished until test R2-3b. 
Finally, R2-4 was intended to assess engine calibration during a full duration test of 159 seconds. 

The main liquid oxidizer orifice and gas generator nozzle orifice were not changed during 
the R2 series. All other orifices were changed after tests R2-2 and R2-3b. Therefore, tests R2-1 
and R2-2 contained identical orifice configurations, as did tests R2-3a and R2-3b. The Engine 5 
calibration baseline was established in test R3-la. Data from this test was used to guide 
selection of an appropriate orifice configuration for use with the high area ratio nozzle. Test 
R3-2b used the same orifice configuration as R2-4 and was used to evaluate Engine 5 calibration 
in response to an orifice change. A summary of engine configuration and test definition 
information for the R2 and R3 test series is provided in Appendix G, Table Gl. 

4.3 Parent Performance Model 

A ROCETS [14] performance model of the MC-1 engine was developed at MSFC to 
support engine design and testing. Fluid conditions in all of the major engine components and 
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flow ducts were modeled using one-dimensional flow physics and empiricisms. This model was 
the nonlinear simulation platform used to generate engine system influence matrices for GDR 
analyses. The ROCETS/MC-1 mainstage model was also modified to run in data reduction (DR) 
mode. This fully nonlinear reduction model, herein designated ROCETS DR, provided the 
conventional data reduction results used to assess the accuracy potential of GDR predictions. 

For computational testing of the GDR strategy, 25 engine system measurements, 
corresponding to performance model variables, were recorded at sampling rates from 25 Hz to 
250 Hz. These measurements and their associated variable names in the ROCETS/MC-1 engine 
model are provided in Appendix G, Table G2. They include four inlet condition measurements 
and 21 internal engine measurements. Fourteen pressures, seven temperatures, two flows, one 
turbopump shaft speed, and the engine thrust comprise the measurement list. 

A total of 22 candidate hardware parameters was originally selected for reduction analysis 
consideration. These are also listed in Table G2 and include four duct/line resistances, four 
valve area ratios, four injector discharge coefficients, two pump head coefficient multipliers, two 
pump torque coefficient multipliers, and one each turbine efficiency multiplier, gas generator 
(GG) exhaust duct orifice coefficient, turbopump friction factor, nozzle discharge coefficient, 
main combustion chamber (MCC) C* efficiency multiplier, and GG liquid oxygen (LOX) inlet 
duct heat transfer rate. 

Only 17 of the 22 candidate hardware parameters were employed in standard data 
reduction analyses performed with the ROCETS DR model. The five parameters not used in 
conventional reductions included the two pump torque multipliers, two main flow valve area 
ratios, and one turbine efficiency. In addition, only 17 of the 21 available internal measurements 
were targeted for matching in conventional reduction analyses. The four measurements not used 
in standard data reduction analyses included the gas generator oxidizer valve (GGOV) inlet 
pressure, the LOX dome temperature, the fuel manifold temperature, and the turbine discharge 
temperature. Hardware parameters and measurements not used in standard data reduction 
analyses are identified in Table G2. Subsystem schematics depicting the approximate flow path 
location of all measurements relative to major hardware components are presented in Appendix 
A, Figures A2 through A5. 

5.0 GDR Pre-anaivsis 

In order to evaluate both the capabilities and limitations of the basic GDR strategy, it is 
necessary to understand the analysis set-up procedure and characteristics of the reduction study 
test data. The next two subsections provide the following information: an outline of the GDR 
analysis setup procedure and a summary of characteristics of the R2 and R3 test series data used 
in reduction analyses. 

5.1 GDR Preprocessing Tasks 

Once candidate measurements and hardware parameters were identified for the MC-1 
engine system, five off-line preprocessing tasks were performed to prepare for a basic GDR 
analysis sequence. These tasks are identified below. 

Preprocessing Tasks 

1. Identify an appropriate baseline state, H 0 C 0 and P 0 
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2. Normalize all measurement variables and hardware parameters 

3. Compute elements of the normalized influence matrices at the base state 

4. Assign weighting factors W for the least squares terms 

5. Assimilate information into appropriately formatted input files. 

Baseline values were selected to closely approximate expected operating conditions of the 
engine system. Baseline values of performance model variables corresponding to internal 
measurements are also required to be an achievable solution of the parent performance model 
with specified baseline values of the hardware parameters and model variables corresponding to 
control/input measurements. The parent performance model for this GDR study was the 
ROCETS/MC-1 mainstage engine model. 

For purposes of this study, baseline values were defined for each MC-1 engine test by the 
following procedure. 

Baseline Definition Procedure 

1.1 The average of data sampled at rates from 25 Hz to 250 Hz (depending on the 
specific measurement) was computed over the time interval from 22 to 23 
seconds for each of the 17 internal measurements and 4 input measurements 
used in standard data reduction. 

1 .2 The averages obtained in the previous step were used as inputs to the ROCETS 
DR model to obtain test specific baseline values for the hardware parameters 

Ho 

1.3 Input variables were assigned the baseline values C 0 indicated below. 


PSVL10 - LOX inlet pressure 47 psia 

TTVL1 0 - LOX inlet temperature 1 64.5 R 
PSRPFV - Fuel inlet pressure 43 psia 

TTRPFV - Fuel inlet temperature 505 R 


1.4 Baseline values P 0 for variables associated with the internal measurements were 
then obtained by running the ROCETS/MC- 1 performance model with 
hardware parameters assigned baseline values H 0 from step 2 and input 
variables assigned the baseline values C 0 indicated in step 3. 

This four-step procedure was adopted to provide a baseline state indicative of the best 
available information regarding engine operation with standard inputs. It was also meant to 
approximate an operational procedure that uses information available from engine acceptance 
tests. Baseline information derived using the above procedure for each MC-1 test is presented in 
Appendix G, Tables G3 and G4. 

The second preprocessing procedure involving normalization of variables is essential for 
robust reduction analysis. The normalization procedure generally affects the parameter selection 
process, the reduction solution H s , and the confidence bounds associated with the solution. The 
preferred scaling procedure involves uncertainty measure normalization. In the absence of 
reliable uncertainty estimates, one approach is to normalize the shift from baseline value for each 
variable by the baseline value for that variable. This is equivalent to uncertainty normalization 
under the assumption that the uncertainty as a percentage of baseline value is the same for all 
variables. Variable and Jacobian normalizations using this assumption are specified below. 
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Normalization Groups 
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The mathematical forms of both the normalized selection process and the normalized 
reduction process in the basic GDR strategy are unchanged using the above groups despite the 
fact that both the selection and solution results may be significantly altered by the rescaling. 

A central difference scheme utilizing ROCETS/MC-1 engine simulation results was 
employed to approximate all Jacobian components. The standard mathematical form of this 
approximation is indicated below. 


P^+SHQ-PiCHoi-SH,-) 

JHo-ij * 2 *SHj 1 J 

In this expression, Pj(H 0 j±5Hj) is the ROCETS/MC-1 output value for dependent variable Pi 
corresponding to an internal measurement with all hardware parameters and input conditions at 
base state conditions except Hj. Baseline incremented values of the hardware parameter, Hj = 
H 0 j±8Hj, are used in the simulation runs that determine values of P, for the numerator difference. 
A similar central differencing procedure was used to obtain the control Jacobian elements Jc 0 -i k - 

The diagonal weighting matrix W provides a means of incorporating information that 
augments the physical and empirical basis of the parent model. If reliable hardware parameter 
uncertainty estimates are available, the diagonal components of the weighting matrix afford an 
effective means of incorporating this information. In this study, the weighting matrix was set 
equal to the identity matrix for all basic GDR analyses. This is consistent with the parameter 
normalizations described in equation (5) and effectively assumes the same level of confidence 
for estimates of all hardware parameters. 

The GDR input file accumulates information required for a reduction run including: 
parameter names and baseline values, parameter normalizing factors, weighting matrix 
components, normalized hardware and control influences, and other parameter designations that 
identity the dimension of the reduction problem. In this study, test data was embedded in the 
GDR input file and processed sequentially, one time slice at a time, to simulate an actual data 
input stream. For real time processing of test or flight data, preloading of data within the input 
file would be replaced by an appropriate interface to the data stream. 

5.2 MC-1 Engine Test Data 

MC-1 engine R2 and R3 test series data for each of the 25 measured quantities identified in 
Table G2 were used. Data gathered during the first five seconds of operation after engine 
ignition were not considered in this study in order to reduce the influence of start up transients. 
One second averaged data over the operating interval from 6 to 23 seconds were used for all tests 
except R2-4. Test R2-4 was a full duration test of 159 seconds, and one second averaged data 
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over the operating interval from 6 to 158 seconds were used for this test. Data for each measured 
variable over the operating interval from 6 to 23 seconds for all tests are displayed in Appendix 
B, Figures B1 through B25. One second averaged data for all seven tests are provided in 
Appendix H, Tables HI through H7. 

General information describing characteristics of the R2 series test data are presented in 
Appendix H, Table H8. For a specific test, the characteristic increase (or decrease) in a 
measured parameter that occurred during the 17 second test interval from 6 to 23 seconds was 
computed by determining the slope of the best fit line relating the measured parameter to time; 
multiplying this result by the 17 second interval length; and then dividing by the test baseline 
value of the parameter. The result is a characteristic increase expressed as a percentage of the 
baseline value. The five-test average increase (or decrease) computed in this fashion is presented 
in Table H8. These values indicate the magnitude and direction of data trends over the 1 7 
second interval from 6 to 23 seconds. The absolute values of these trends generally fell in the 
interval from one to three percent. A significant exception was the inlet LOX pressure which 
exhibited an average increase of more than 6%. 

As mentioned earlier, the same orifice set - number 5 - was used in tests R2-1 and R2-2. For 
each measured parameter, the absolute value of the difference between the test R2-1 observed 
value at a given time slice, and the test R2-2 observed value at the same time slice, averaged over 
all time slices, provides an indication of the normal dispersion of the measured parameter for a 
constant configuration engine system. Similar dispersion information for a constant 
configuration engine system using orifice set 6 was obtained by examining data from tests R2-3a 
and R2-3b. These test-to-test dispersion measures are presented in Table H8 for each measured 
parameter. The observed differences are caused by inlet condition, hardware function, and 
measurement system variation. For example, the effects of engine fuel inlet pressure (PSRPFV) 
variation are clearly mirrored by the fuel pump inlet pressure (PSVLOO) difference in the orifice 
set 5 test sequence R2-1 R2-2. The large difference in turbine discharge temperature 
(TTHTGD) for the orifice set 5 test sequence was caused by a temperature sensor failure during 
testR2-l. 

The effects of inlet condition variation can be removed by adjusting values of all measured 
internal parameters for standard conditions, i.e., expected values given baseline inlet conditions 
and baseline hardware conditions. Using the linearized performance model relations from 
equation (1), expected values of the measured internal parameters adjusted for standard baseline 
inlet conditions can be estimated by the relation below. 

Pstd = P meas ' Jco (C m eas — C 0 ) (8) 

Plots of R2 series test data adjusted for standard inlet conditions are shown in Appendix C, 
Figures Cl through C21. Examination of these plots indicates a significant reduction in the rates 
of increase (decrease) of most measured internal parameters over the sampled range from 6 to 23 
seconds. These rate reductions are confirmed by the average R2 series rise (drop) values for 
standard inputs shown in Appendix H, Table H9. This suggests that much of the time dependent 
variation in measured parameters was due to inlet condition variation over the course of the test. 

Orifice set 5 average absolute differences for the test sequence R2-1 R2-2 were reduced for 
all internal measurements adjusted for standard inlet conditions as shown in Table H9. However, 
orifice set 6 differences, associated with the test sequence R2-3a R2-3b, were generally greater 
after correction to a standard input state. For each internal measurement parameter, the observed 
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test-to-test average absolute difference after adjustment for standard inlet conditions was caused 
primarily by hardware function and measurement system variation. The oxygen bleed valve leak 
that occurred in test R2-3a represents an additional unknown boundary condition variation that 
impacted orifice set 6 differences. 



In order to facilitate comparison with accepted reduction results obtained from the parent 
ROCETS DR model, the same 17 hardware parameters used in the parent reduction process were 
used in GDR analyses of the MC-1 engine. These 17 parameters have been identified previously 
and correspond to the unshaded hardware parameters shown in Appendix G, Table G2. 

After preprocessing operations were completed, the automated subset selection process was 
performed to eliminate redundant or ill-scaled measurements, thereby reducing the L 2 condition 
number of the hardware Jacobian matrix to an acceptable value. In all cases the first four 
measurements eliminated by the automated selection process were as indicated below. 


TTVL05 Fuel manifold temperature 

TTVL14 LOX dome temperature 

TTHTGD Turbine discharge temperature 

PSVL15 Gas generator oxidizer valve (GGOV) inlet pressure 

These were precisely the four measurements eliminated by expert selection prior to ROCETS 
DR analyses. After elimination of these four measurement variables, the resulting Jacobian 
condition number was not small enough to guarantee stability. However, basic GDR analyses 
were run using the 1 7 remaining internal measurements and 4 inlet measurements to facilitate 
direct comparison to the accepted ROCETS DR standard. 

Typical GDR and ROCETS data reduction results for test R2-1 are displayed in Appendix 
D, Figures D1-D17. These figures contain both GDR and ROCETS time history data reduction 
results for each of the 17 selected hardware parameters. Differences between GDR and 
ROCETS DR results are generally observed to be small over the entire test time interval. The 
largest differences are associated with the fuel pump inlet line resistance RKFL1 and the gas 
generator fuel valve area ration XMGGKO. 

The difference between GDR and ROCETS DR predictions was computed for each 
hardware parameter at each one second time slice. Absolute values of these differences were 
averaged over all time slices in the test interval for each hardware parameter. This information is 
displayed in Appendix I, Table II . Average absolute deviations greater than 0.5% of baseline 
were confined to test R2-1 and test R2-2. This indicates that the basic GDR procedure was 
capable of providing a consistent approximation of hardware behavior that agrees well with the 
recognized standard for MC-1 engine data reduction. However, average absolute deviations in 
excess of one percent for two parameters in each of tests R2-1 and R2-1 plus the large Jacobian 
condition number were concerns. 

It should be noted that the primary motivation for constructing a simplified representation of 
the performance model is to improve analysis speed to support health monitoring. The basic 
GDR results for the MC-1 engine were processed at a rate greater than 200 time slice reductions 
per second on an 800 MHz Pentium III processor. This qualifies the procedure as a real-time 
analysis tool. There are other advantages to the simplified GDR representation of the engine 
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related to stability and accuracy characterization that are addressed later in this report. The next 
section documents enhancements considered for the basic GDR procedure, which will henceforth 
be identified as GDRA, in order to improve computational speed, accuracy, and stability. 

7.0 GDR Modifications 

Three types of modifications to the basic GDR strategy were explored in efforts to improve 
the effective solution range of the method as well as the computational speed, accuracy, and 
stability of the solution process. 

Types of GDR modifications 

1. Modifications of the subset selection scheme 

2. Modifications of the solution procedure 

3. Modifications of the problem formulation 

Specific changes investigated during the course of this effort are described in the following 
subsections. 

7.1 Solution Procedure Using Singular Value Decomposition 

A very efficient procedure for determining the minimum 2-norm solution to a linear system 
is based on the singular value decomposition (SVD) of the coefficient matrix. This procedure, 
applied to the balance relations described in equation (1), requires computation of the SVD of the 
hardware Jacobian matrix. The form of the decomposition is described below [see, e.g. 6]. 

U T J Ho V = I (9) 

where 

U is a mxm orthogonal matrix with column partitioning U = [ uj . . . u m ] 

V is a nxn orthogonal matrix with column partitioning V = [ Vi . . . v„ ] 

£ = diagOi . . . ct p ) is the diagonal matrix of singular values Oj arranged in 
descending order of magnitude 
p = min(m,n) 

Using information from the decomposition, the minimum 2 norm solution to the system 

JhoAH = b = AP-JcoAC (10) 


can be computed directly using the relation 


AH = 



i=l 


( 11 ) 


where r < p is the estimated rank of Jh 0 - 

Equation (1 1) is obtained using a generalized inverse of the hardware influence matrix to 
isolate the hardware solution. Once the SVD of the Jacobian Jh 0 is computed, hardware 
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solutions may be obtained very quickly from the direct computation of AH prescribed by 
equation (1 1). A version of GDR based on the SVD decomposition procedure described in 
equations (9) through (1 1) was implemented and will be referred to as GDRB. 

Computational experiments were performed with GDRB to determine efficiency and to 
assess potential for reducing subset selection overhead. For near rank deficient Jacobians, 
reduction in the rank estimate r, which corresponds to elimination of near singular directions, can 
be used instead of parameter elimination to reduce redundancy and scaling problems. This 
effectively reduces the need for a separate subset selection process. Unfortunately, rank 
reduction has two significant drawbacks that disqualified it as a subset selection substitute. 

1. If a rank estimate r < p is used in equation (1 1), hardware shifts are no longer a 
continuous function of the data 

2. Realistic single parameter shifts may not be recoverable upon elimination of singular 
direction(s) implicit with rank reduction. 

Because of the computational overhead added to complete the initial SVD, improvements in 
GDRB reduction cycle speed were negligible for short duration test runs. The average reduction 
cycle rate for test R2-4 - the full duration test of 159 seconds - was increased by approximately 
30%. Coupled with the QR-Pi subset selection routine, GDRB reduction solutions were 
essentially the same as those returned by the basic GDR procedure. 

7.2 Partial Second Order Formulation and Solution Procedure 

The primary range and accuracy restrictions of the GDR procedure are imposed by the 
model linearization leading to equation (1). There are many paths available to relax the linear 
approximation. Perhaps the most direct approach is to incorporate higher order effects in a 
sequence beginning with second order terms. For a system with 17 defining hardware 
parameters, there are 153 distinct second order terms — 136 covanant terms and 17 second order 
terms of the form (AHj) 2 . If component hardware performance parameters are defined such that 
there is little interaction between the parameters, i.e., the performance parameters are 
independent, then the covariant terms are negligible. 

Although the hardware parameters used for MC-1 engine reduction analyses were not 
strictly independent, a second order extension of GDR adding only terms of the form (AHj) was 
implemented and tested. This extension, designated GDRC, used the approximate second order 
balance relations below to define the reduction solution. 


AP « J H o AH + J Co AC + \ D Ho diag(AH) AH + \ D Co diag(AC) AC 


Dho - 


A five point central difference scheme was employed to compute elements of the second 
order influence matrices Dho and I)( 0 - The measurement selection process was assumed static — 
utilizing only the first order hardware Jacobian at the base state in the QR-Pi factorization that 
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identified measurements to eliminate. A Newton type iterative method [15 ] was used to solve 
balance relations (12) reorganized in the form below. 


Jho + 2 D Ho diag(AH) 


AH = AP - Jco AC - j D Co diag(AC) AC 


( 14 ) 


Because deviations from standard DR results returned by GDRA were larger for tests R2-1 
and R2-2 (see Appendix I, Table II), nonlinear GDRC analyses were performed using only R2-1 
and R2-2 test data. Comparisons to GDRA results are given in Appendix I, Table 12. 

It is obvious from Table 12 that predicted deviations from standard DR results were not 
consistently improved using the nonlinear GDRC formulation. In addition, the deviation of 
GDRC predictions from standard DR results for the GG fuel valve area multiplier XMGGKO 
were significantly greater than those returned by GDRA. Further testing was not performed 
given the greater computational overhead and lack of significant improvement in solution 
accuracy. 


7.3 Formulation with Hardware Parameter Functions 

Modification of the original hardware parameter set was indicated based on the following 
observations. 


1 . Influences of the main fuel line calibrating resistance RCALMF were large 
relative to hardware specific resistances along the main fuel line. 

2. Influences of the main oxidizer line calibrating resistance RCALMO were large 
relative to hardware specific resistances along the main oxidizer line. 

3. Shaft friction factor FRICFACT shifts were very large as a percent of baseline 
at the reduction solution. This indicated that FRICFACT was not well scaled 
and could artificially dominate a linear systems solution process. 

4. All valve area multipliers and injector discharge coefficients were proportional 
to one over a driving pressure difference squared (i.e., valve XM’s and injector 
CD’s oc 1/AP 2 ). Therefore, pressure dependence on XM’s and Cd’s is 
inherently nonlinear. 

These observations motivated a revision of the set of hardware parameters used for MC-1 
engine reduction analyses as indicated in Table 1 below. The complete set of hardware 
parameters used in data reduction analyses after incorporating the revisions below is 
given in Appendix G, Table G5. 
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Table 1 MC-1 engine revised hardware parameters 


Original Hdwe 
parameter 

Modified hdwe 
parameter 

Modified hdwe description 

Spanning Pressures 

CDGGKI 

RGGKI 

GG fuel injector equivalent resistance 

PSVL09 - PTHTGI 

CDGGOI 

RGGOI 

GG oxidizer injector equivalent resistance 

PSVL18 - PTHTGI 

CDKINJ 

RKINJ 

MCC fuel injector equivalent resistance 

PTVL05 - PTMCHY 

CDOINJ 

ROINJ 

MCC LOX injector equivalent resistance 

PTVL14 - PTMCHY 

XMMCKO and 




RCALMF 

RMMCRP 

Main fuel line equivalent resistance 

PSVL01 - PTVL05 

XMMCOO and 




RCALMO 

RMMCOX 

Main LOX line equivalent resistance 

PSOXDS - PSVL13 

XMGGKO 

RMGGRP 

GG fuel line equivalent resistance 

PSVL01 - PTVL09 

XMGGOO 

RMGGOX 

GG LOX line equivalent resistance 

PSVL15 - PSVL18 

FRICFACT 

PWRFACT 

Turbopump power factor 



PWRFACT in the above list can be thought of as the ratio of actual pump shaft delivered 
power to actual turbine shaft supplied power. It combines the effects of turbine efficiency 
multiplier, pump torque multipliers, and shaft friction factor such that the following 
proportionality applies 


PWRFACT oc 


1-FRICFACT 

ETAMHTGT * (2 pump equivalent torque multiplier) 


(15) 


The switch to equivalent resistances is a way of converting the original hardware parameter 
set to a more natural parameter set composed of linearly independent nonlinear functions of the 
original parameters. Since resistance is by definition proportional to spanning pressure drop, 
measured pressures are more naturally represented as linear functions of the equivalent 
resistances. 

In general terms, the data reduction problem with the modified hardware parameters can 
be represented as follows. 


Find Af(H) such that 

AP ~ Jf(Ho) Af(H) + Jco AC 

where 

f(H) = [ f,(H) . . . f n (H) ] T Af(H) = [ f,(H)-f,(H 0 ) . . . f„(H)-f„(H 0 ) ] T 


Jf(Ho) - 



For the MC-1 engine, the functions f(H) represent the modified hardware parameter set 
identified in Table G5. 


(16) 

(17) 

(18) 


7.3.1 GDRA and ROCETS DR Comparisons 
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GDRA analyses of the MC-1 engine were performed using the revised hardware parameter 
set listed in Table G5. The deviation of GDR predicted values from standard ROCETS DR 
results, converted to display values of the revised hardware parameters, is provided in Appendix 
I, Table 13. The agreement is observed to be excellent. Over the seven test R2 and R3 series, 
only the inlet fuel resistance RKFL1 in test R2-1 exhibited an average deviation greater than one 
percent. Time slice reduction cycles were performed at a rate greater than 200 per second on a 
system equipped with an 800 MHz Pentium III processor. The accuracy of these reduction 
results coupled with high cycle speed qualifies the process as a real time monitoring candidate 
for nominal MC-1 engine operation. 

7.3.2 Single Source Anomaly Resolution Capability 

In order to assess GDR capability to identify single source anomalies, a series of simulations 
were performed with single source hardware anomaly input. Simulated test data for each single 
source anomaly was obtained as output of ROCETS/MC-1 runs with one hardware parameter 
modified from its baseline value. This data was provided for GDRA reduction analyses to 
determine if the anomaly would be identified. Results are presented in Appendix E, Figures El 
through E20. Each figure identifies the magnitude of the input anomaly and the distribution of 
causation allocated by each of ROCETS data reduction and GDRA using simulated test data. 
Two GDRA results are presented, one that employed the standard 17 measurement set also used 
to derive ROCETS DR results, and one using only 16 measurements as an example of lost 
measurement effect on anomaly resolution capability. The additional measurement eliminated 
for the 16 measurement case was the LOX GG inlet temperature TTVL18. 

High fidelity anomaly recovery is generally observed. The following comments address 
specific characteristics of the simulation results. 

1 . Figure E 1 3 . A fuel pump torque multiplier (TRQMKPMP) anomaly reduction is 
recovered as a clearly identifiable power factor (PWRFACT) increase by all reduction 
procedures. This is consistent with the definition of power factor. 

2. Figure El 6. Same comments as number 1 for a turbine efficiency multiplier 
(ETAMHTGT) reduction. 

3. Figure E20. LOX inlet duct heat transfer (QDOTVL18) shift was not recovered by 
GDRA using 16 measurements. This is because the GG LOX inlet duct temperature 
TTVL1 8 is a direct indicator of duct heat transfer. When this measurement is lost, the 
ability to identify heat transfer rate is also lost. This is an example of a very specific loss 
of discrimination effect. Other effects of TTVL18 measurement loss on anomaly 
resolution capability are negligible. 

4. Figure El 4. An oxygen pump torque multiplier (TRQMOPMP) reduction was allocated 
as various effects: a main oxygen line equivalent resistance (RMMCOX) increase, a gas 
generator oxygen line equivalent resistance increase (RMGGOX), and a power factor 
(PWRFACT) increase by all reduction methods, and a LOX duct heat transfer 
(QDOTVL18) reduction by both ROCETS DR and GDR using 17 measurements. This 
was the only case in which no reduction procedure, including the parent ROCETS DR 
procedure, could allocate a single source anomaly to a single predominant reduction 
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parameter. This indicates a need to modify the measurement suite and/or the parent 
performance model to allow clear identification of anomaly operation. 

It should be noted that shifts in both the fuel pump torque multiplier (TRQMKPMP) and the 
turbine efficiency multiplier (ETAMHTGT) have a clear power factor (PWRFACT) signature. 
This suggests that shifts in these parameters can be detected but not discriminated using the 
revised reduction parameter set. The oxygen pump torque multiplier (TRQMOPMP) is the only 
parameter that did not have a clear anomaly signature. 

7.3.3 Sensor Failure Response 

A study was performed to determine the effects of single sensor failures on reduction 
analysis predictions using the revised set of hardware parameters. GDRA solutions with 
individual sensors eliminated were obtained and compared to standard ROCETS DR results 
using all 17 measurements. Results using test R2-4 data are reported in Appendix I, Table 14, in 
the form of average absolute deviations from ROCETS DR results as a percentage of baseline 
hardware parameter value. 

In each column of Table 14, measurements eliminated by subset selection are identified in 
addition to the measurement eliminated to simulate sensor failure. Hardware parameters with 
deviations in excess of 2% are identified as unreliable with the indicated sensor failure. 
Parameters with deviations between 1% and 2% are identified as only marginally reliable with 
the indicated sensor failure. Results indicate that the total LOX flow (WOXTOTL) and the 
turbine inlet temperature (TTHTGI) are the most critical for reduction accuracy. 

7.3.4 Results with Flight Measurement Suite 

Only five internal measurements in addition to the four inlet measurements were planned 
for MC-1 flight engines. The flight measurement set is identified in Appendix G, Table G6, and 
includes one temperature at the turbine inlet (TTHTGI), and four pressures: two pump discharge 
pressures (PSOXDS and PSVL01), the gas generator pressure (PTHTGI), and the main chamber 
pressure (PTMCHY). No flow measurements were anticipated for the flight engine. 

The candidate hardware parameter set was reduced to 12 by combining resistances between 
pressure measurements and by eliminating the nozzle discharge coefficient characterized 
primarily by thrust. The set of candidate hardware parameters for flight data reductions is 
defined in Table G6. 

Three of the flight candidate hardware parameters were further eliminated because relevant 
temperature measures were not available. These included the two pump torque multipliers 
(TRQMKPMP and TRQMOPMP) and the turbine efficiency multiplier (ETAMHTGT). GDRA 
analyses were performed using data from test R2-4 for five internal measurements and four inlet 
condition measurements to estimate the nine remaining hardware parameters. Results of these 
analyses are compared to ROCETS DR results with a full suite of 1 7 measurements in Appendix 
F, Figures 1 through 9. Summary observations are provided below. 

1. Loss of detail discrimination is evident for main fuel line series resistance R3MCRP in 
Figure FI. 
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2. Some loss of detail and trending is observed for main oxidizer line series resistance 
R3MCOX in Figure F2. 

3. Little or no resolution loss is observed for GG inlet series resistances R3GGRP in Figure 
F3 and R3GGOX in Figure F4. 

4. Little discrimination loss and only slight tail end trending loss is observed for the fuel 
pump head coefficient multiplier in Figure F5. 

5. Some trending loss is observed for the LOX pump head coefficient multiplier in Figure 
F6. 

6. Very slight trending loss is observed for the power factor term PWRFACT in Figure F7 . 

7. Some trending loss is observed for the MCC efficiency multiplier in Figure F8. 

8. Heat transfer effects QDOTVL1 8 were not recovered at all as indicated in Figure F9. 

Average and maximum absolute deviations of GDRA results, obtained using only flight 
measurements, from standard ROCETS DR predictions, obtained using a full suite of 17 
measurements, are provided in Appendix I, Table 15. Considering the limited measurement 
input, results were judged to be quite good with only one parameter (QDOTVL18) deviating 
over 2 % from its standard value. 

In an attempt to correct some of the detail and trending losses apparent in Figures FI 
through F9, fixed LOX and fuel flow rates obtained by averaging flow data over the time interval 
from 20 to 23 seconds were added to the flight data for reduction analysis. This was intended to 
simulate the use of best available flow information from acceptance testing. Reduction results 
obtained using the fixed flow rates are presented in Appendix F, Figures F10 through FI 8. 
Summary observations are provided below. 

1 . Most of the discrimination loss reported for R3MCRP was recovered as shown in Figure 
F10. 

2. Most of the trending loss reported for R3MCOX was recovered as shown in Figure FI 1 . 

3. Most of the trending loss reported for PSIMOPMP was recovered as shown in Figure 
FI 5. 

4. Trending losses were exacerbated for PWRFACT as shown in Figure FI 6. 

5. Most of the trending loss reported for ECSMMCHB was recovered as shown in Figure 
FI 7. Detail discrimination was reduced however. 

6. No improvement in QDOTVL18 recovery was observed as shown in Figure FI 8. 

Average and maximum absolute deviations of GDRA results, obtained using flight 
measurements with 20-23 second average flows, from standard ROCETS DR predictions, 
obtained using a full suite of 17 measurements, are also given in Appendix I, Table 15. With the 
exception of QDOTVL18, some improvement is evident especially in the standard deviation of 
the absolute differences. GDRA reduction result using flight measurements augmented with 
fixed average flows indicative of actual engine operation, such as available from acceptance 
tests, were evaluated as reasonable for real-time monitoring applications. 

7.4 Subset Selection Procedures 

One of the key components of the GDR strategy is a parameter selection strategy. 

Automated parameter subset selection is applied primarily to eliminate measurements and 


18 


stabilize ill-conditioned underdetermined systems. In the basic GDR method, this is 
accomplished using a QR-Pi factorization of the transpose of the hardware Jacobian matrix. 
Although heuristic, this method appears to effectively isolate and eliminate measurement 
redundancy and scaling problems. 

Automated hardware parameter selection is inherently more difficult. For robust engine 
health monitoring, hardware parameters should be selected based primarily on 
failure/degradation risk and impact. Failure/degradation potential can only be established by 
examining the experience base of like components and engine systems. Lacking effective 
estimates, hardware parameter elimination is difficult to justify and can be detrimental if 
components with significant failure risk are eliminated. 

Despite the above observations, two automated procedures for selecting both hardware and 
measurement variable eliminations were constructed and implemented. One was based on 
sequential QR-Pi factorization of the hardware Jacobian and its transpose. The second was 
based on sequential condition number reduction. In both cases, the order of hardware and 
measurement eliminations is user specified. 

Computational experience with the dual selection strategies is still limited. Typical results, 
obtained using the QR-Pi dual elimination algorithm to select 5 measurements from the 
available suite of measurements for the MC-1 engine, are presented in Appendix I, Table 16. 

The leftmost shaded section of the table displays the flight internal measurement set specified for 
the MC-1 engine. The remainder of the table presents results of four automated subset selection 
sequences differing by user specified front-end eliminations and measurement-hardware 
elimination sequence. Candidate measurement variables, before any elimination, were the 21 
MC-1 internal measurements identified in Appendix G, Table G2. Candidate hardware 
parameters, before any elimination, included only members of the revised hardware set defined 
in Table G5. 

The procedure designated “Auto-Selection 1” in Table 16 was a measurement only 
elimination sequence beginning with four user defined hardware eliminations and one 
measurement parameter - engine thrust - elimination. The automated procedure eliminated 1 5 
additional measurements. The five remaining measurements variables were the LOX pump 
discharge pressure (PSOXDS), the turbine inlet temperature (TTHTGI), both the fuel and LOX 
total flows (WRPTOTL and WOXTOTL), and the turbine discharge pressure (PTVL22). Only 
two of the five selected measurements is included in the flight list. There are many possible 
reasons for this difference that have nothing to do with the mathematical characteristics of the 
hardware Jacobian. These include a host of practical consideration including cost, ease of 
installation and maintenance, reliability, life, uncertainty, etc. 

“Auto-Selection 2” is another measurement only elimination sequence, but beginning with 
three user defined eliminations including both propellant flow rates. The selection process 
identified five measurements from a list that does not include the propellant flows. The five 
selected measurements include three on the flight list. 

“Auto-Selection 3” is a repeated measurement-hardware elimination sequence that selects 
five measurements and five compatible hardware parameters. The selected measurements are 
the same as identified in auto-selection 2. The selected hardware parameters include critical 
turbopump (PSIMKPMP, PSIMOPMP, and PWRFACT), chamber (ECSMMCHB), and gas 
generator nozzle (CDGGNZ) characteristics. 

“Auto-selection 4” is another repeated measurement-hardware elimination sequence that, 
because of user eliminations, begins with a smaller set of hardware parameters. The parameters 
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selections are similar to those of “Auto-Selection 3” and the reader is referred to Table 16 for 
detailed information. 

It is important to understand that sequential selection procedures do not identify globally 
optimum subsets in general. This is a major drawback and suggests the need for computational 
testing to validate the suitability of parameter selections. 

The mathematical and statistical literature dealing with subset selection is extensive [see 
e.g., 6, 15, 16, 17, 18]. A robust selection strategy must incorporate risk and uncertainty 
information that goes beyond classical analytical considerations of stability. Uncertainty scaling 
is one natural method of incorporating uncertainty information within the selection procedure; 
however, the assignment of risk and uncertainty estimates appropriate for MC-1 engine operation 
was beyond the scope of this effort. 

8.0 Summary of Research Effort 

The following procedures were developed and implemented during the course of this 
research effort. 

1 . An efficient generalized inverse solution procedure based on singular value 
decomposition of the hardware Jacobian. This solution procedure was incorporated in 
the code version identified as GDRB. 

2. A version of GDR that incorporates a restricted second order approximation model of 
engine performance and a Newton type solution scheme. This approximation model and 
solution procedure was incorporated in the code version designated GDRC. 

3. A redefined hardware parameter formulation for the MC-1 engine that improves the 
linear system approximation and a general formulation that characterizes hardware 
combination schemes 

4. Two dual sequential subset selection schemes that identify computationally compatible 
combinations of hardware parameters and measurement variables for data reduction. 

The following GDR performance measures for MC-1 engine system reduction analysis were 
obtained. 

1 . GDRA reduction analysis cycle rate >200 cycles per second, GDRB > 250 cycles per 
second, GDRC >50 cycles per second 

2. Average hardware parameter deviation from full nonlinear parent model data reduction 
results using full measurement suite and revised hardware list, < 0.5% except for isolated 
cases (Table 13) 

3. Single source anomaly recapture capability with full measurement suite and revised 
hardware list (based on simulations), >90% of parent model signature with one possible 
exception (Figures El through E20) 

4. Qualified integrity of reduction hardware parameter recovery for single sensor losses and 
identified loss of discrimination parameters for each sensor loss (Table 14) 

5. Average deviation of GDRA hardware predictions, using only the flight measurement set, 
from full nonlinear parent model DR results using full measurement suite, <1% with two 
exceptions, <2% with one exception (Table 15) 
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9.0 Assessments and Recommendations 


The following assessments are specific for the MC-1 engine system instrumented with the 
standard measurement suite (Table G2) and GDRA using the revised hardware parameter set 
(Table G5). 

1. Assuming sensor integrity, GDR would be an effective real time monitoring tool for the 
MC-1 engine system. 

2. GDR is capable of detecting and isolating most single source hardware anomalies and 
should be capable of resolving many multiple source anomalies. GDR cannot identify 
oxygen pump efficiency anomalies and cannot discriminate between efficiency losses in 
the turbine or fuel pump. 

3. GDR can reliably track the behavior of most hardware parameters in the presence of a 
single sensor failure if the failure is recognized. 

4. GDR can predict nominal trends of flight revised hardware parameters (Table G6) given 
the five specific internal flight measurements identified in Table G6 and good estimates 
of engine component propellant flows. Reliable anomaly resolution capability using only 
the designated flight measurements is unlikely. 

The following are general assessments and recommendations based on results of this study. 

1 . The QR-Pi dual subset selection routine developed during this effort provides an 
effective sequential method of specifying compatible measurement and hardware sets for 
reduction analyses. It is limited to considerations of scale and redundancy based on 
information from the best available system model. 

2. Significant improvement in selection of appropriate measurements for robust health 
monitoring will require reliable estimates of both measurement and model uncertainty. 

3. If reasonable uncertainty estimates of measurement and model uncertainty are available, 
GDR should be modified to incorporate uncertainty scaling. This would improve both 
the subset selection process and the solution procedure. 

4. A globally optimum measurement and hardware selection procedure should be 
developed. 

5. The GDR development path is incremental. Reliable extensions will require application 
to an engine system with a larger controlled operating range. 

6. GDR predictions are susceptible to data bias. GDR should be used in conjunction with a 
reliable sensor qualification algorithm, or a qualification procedure should be 
incorporated in GDR. 

7. A simulation platform to study the effects of noise and measurement uncertainty on 
prediction reliability should be developed. 
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Appendix A 


MC-1 engine propellant flow schematics 


Figure A1 MC-1 Engine Conceptual Schematic 





Figure A2 MC-1 Engine Oxygen System 
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Figure A3 MC-1 Engine Fuel System 






Figure A4 MC-1 Engine GG/Turbine System 
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Figure A5 MC-1 Engine MCC/Nozzle System 




Appendix B 


MC-1 engine 

Temporal plots of R2 and R3 test series 
One second averaged data 
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Figure B2 TTVL10 one second average test data 
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Figure B3 PSRPFV one second average test data 
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Figure B5 PSOXDS one second average test data 
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Figure B7 PTVL14 one second average test data 
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Figure B8 PSVL15 one second average test data 
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Figure B9 PTVL18 one second average test data 
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Figure BIO TTVL14 one second average test data 
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Figure B1 1 TTVL18 one second average test data 
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Figure B12 WOXTOTL one second average test data 
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Figure B13 PSVLOO one second average test data 
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Figure B14 PSVL01 one second average test data 
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Figure B15 PTVL05 one second average test data 
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Figure B16 PTVL09 one second average test data 
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Figure B17 TTVL05 one second average test data 
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Figure B18 WRPTOTL one second average test data 
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Figure B19 PTHTGI one second average test data 
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Figure B20 PTVL22 one second average test data 
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Figure B21 TTHTGI one second average test data 
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Figure B22 TTHTGD one second average test data 
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Figure B23 PTMCHY one second average test data 
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Figure B24 SNSHFT one second average test data 
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Figure B25 FT15A one second average test data 
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Appendix C 


MC-1 engine 

Temporal plots of R2 and R3 test series 
One second averaged data adjusted 
to standard inlet conditions 


Figure Cl PSOXDS test data adjusted to standard inlet conditions 
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Figure C2 PSVL13 test data adjusted to standard inlet conditions 
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Figure C3 PTVL14 test data adjusted to standard inlet conditions 
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Figure C4 PSVL15 test data adjusted to standard inlet conditions 
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Figure C5 PTVL18 test data adjusted to standard inlet conditions 
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Figure C6 TTVL14 test data adjusted to standard inlet conditions 
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Figure C7 TTVL18 test data adjusted to standard inlet conditions 
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Figure C8 WOXTOTL test data adjusted to standard inlet conditions 
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Figure C9 PSVLOO test data adjusted to standard inlet conditions 
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Figure CIO PSVLOI test data adjusted to standard inlet conditions 
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Figure Cl 1 PTVL05 test data adjusted to standard inlet conditions 
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Figure Cl 2 PTVL09 test data adjusted to standard inlet conditions 
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Figure Cl 3 TTVL05 test data adjusted to standard inlet conditions 
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Figure Cl 4 WRPTOTL test data adjusted to standard inlet conditions 
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Figure Cl 5 PTHTGI test data adjusted to standard inlet conditions 
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Figure Cl 6 PTVL22 test data adjusted to standard inlet conditions 
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Figure Cl 7 TTHTGI test data adjusted to standard inlet conditions 
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Figure Cl 8 TTHTGD test data adjusted to standard inlet conditions 
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Figure Cl 9 PTMCHY test data adjusted to standard inlet conditions 
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Figure C20 SNSHFT test data adjusted to standard inlet conditions 
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Figure C21 FT15A test data adjusted to standard inlet conditions 
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Appendix D 


MC-lengine 

Temporal plots for test R2-1 
Hardware parameter variation with the 
standard reduction variable set 


Figure D1 Comparison of GDR and ROCETS/DR results 
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Figure D2 Comparison of GDR and ROCETS/DR results 
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Figure D3 Comparison of GDR and ROCETS/DR results 
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Figure D4 Comparison of GDR and ROCETS/DR results 
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Figure D5 Comparison of GDR and ROCETS/DR results 
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Figure D6 Comparison of GDR and ROCETS/DR results 
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Figure D7 Comparison of GDR and ROCETS/DR results 
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Figure D8 Comparison of GDR and ROCETS/DR results 
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Figure D9 Comparison of GDR and ROCETS/DR results 
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Figure DIO Comparison of GDR and ROCETS/DR results 
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Figure Dll Comparison of GDR and ROCETS/DR results 
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Figure D12 Comparison of GDR and ROCETS/DR results 
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Figure D13 Comparison of GDR and ROCETS/DR results 
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Figure D14 Comparison of GDR and ROCETS/DR results 
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Figure D15 Comparison of GDR and ROCETS/DR results 
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Figure D16 Comparison of GDR and ROCETS/DR results 
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Figure D17 Comparison of GDR and ROCETS/DR results 
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Appendix E 


MC-1 engine 

Allocation of single source anomaly 
effects by data reduction 


Figure El RMMCRP single source anomaly allocation 
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Figure E2 RMMCOX single source anomaly allocation 
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Figure E3 RKFL1 single source anomaly allocation 
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Figure E4 R0LN1 single source anomaly allocation 
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Figure E5 RMGGRP single source anomaly allocation 
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Figure E8 RGGOI single source anomaly allocation 
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Figure E10 ROINJ single source anomaly allocation 
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Figure Ell PSIMKPMP single source anomaly allocation 
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Figure E12 PSIMOPMP single source anomaly allocation 



8HA10QD 

IGVddMd 

X099V\ld 

dd99IAId 

dlAldOlAIDdl 

dlAldXI/MOdl 

miod 

XOOIAIIAId 

daoi^iiAia 

dIAIdOIAIISd 

dlAldXIAIISd 

191HIAIV13 

9H0I/M1AIS03 

IZONOO 

zNOGao 

rNiod 

rNixd 

I099d 

IX99d 


£ 

IE 

co 


E-12 


Reduction Assignment 




8HA10QD 


IQVddMd 

XOOOlAld 

ddooiAia 

dlAJdOlAIDdl 

dl/MdXlAIDdl 

UMlOd 

nd»d 

XOOlAIIAld 

ddOlAllAld 

dIAIdOIAIISd 

dIAIdXIAIISd 

101HIAIV13 

9H0IAIIAIS03 

1ZONQO 


ZNOOQO 


Reduction Assianment 






8HAiOaO 


lovd^Md 

XOOOlAld 

ddOOIAId 

dlAldOIAlOdl 

dlAldXIAIOdl 

1-NlOd 

HdXd 

XOOIAIWd 

ddOlAllAld 

dIAIdOiAJ ISd 

dIAIdXIAIISd 

191HIAIV13 

9H01AI1AIS03 

1ZONQO 


ZN99Q0 


Reduction Assignment 





8HA10aD 


lovdyMd 

xooowu 

ddOOlAld 

dlAldOlAlDdl 

dlAldXIAIOdl 

miod 

ndxd 

XOOlAllAld 

ddOIAIIMd 

dIAIdOIAIISd 

dIAIdXIAIISd 

101H1AIV13 

aHOlAllAISOd 

IZONdO 


ZN90Q0 


Reduction Assianment 





8 nAioao 


lOVddMd 

xoooiAia 

daoowd 

dlAldOIAIOdl 

dlAldXWOdl 

INIOU 

ndxd 

XOOlAllAld 

ddOlAIIMd 

dIAldOlAIISd 

dIAIdXIAIISd 

191HIAIV13 

gHOIAIWSOd 

1ZONQO 


ZNOoao 


Reduction Assianment 





9nA10QD 


lOVddMd 

XOOOlAJd 

dd99IMd 

dl/\ldOI/\IOd± 

dlAldXlAIOdl 

miod 

nd>id 

XOOIAIlAld 

ddOlMlAld 

dlAldOlAIISd 

dlAldXWISd 

JL91HIAIV13 

9H01MNS03 

1ZONQO 


ZN99Q0 


Reduction Assianment 






Figure El 8 CDNOZL single source anomaly allocation 
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Figure E19 ECSMMCHB single source anomaly allocation 
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Appendix F 


MC-1 engine 

Comparison of standard ROCETS DR andGDRA using 

only flight measurements 



Figure FI Reduction results with flight measurement set 
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Figure F2 Reduction results with flight measurement set 
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Figure F3 Reduction results with flight measurement set 

R3GGRP 
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Figure F4 Reduction results with flight measurement set 
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Figure F5 Reduction results with flight measurement set 
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Figure F6 Reduction results with flight measurement set 
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Figure F7 Reduction results with flight measurement set 

PWRFACT 
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Figure F8 Reduction results with flight measurement set 

ECSMMCHB 
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Figure F9 Reduction results with flight measurement set 

QDOTVL18 
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Figure F10 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure F11 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure FI 2 Reduction results with flight measurement set and 20-23 second avg flows 

R3GGRP 
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Figure FI 3 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure FI 4 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure FI 5 Reduction results with flight measurement set and 20-23 second avg flows 

PSIMOPMP 



F-15 


Time (sec) 





Figure FI 6 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure FI 7 Reduction results with flight measurement set and 20-23 second avg flows 
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Figure FI 8 Reduction results with flight measurement set and 20-23 second avg flows 
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Appendix G 


MC-1 engine 
R2 and R3 test series 
Information for data reduction 



Table G1 Description of MC-1 engine test sequence used as study data source 






G-2 



Table G3 Baseline values for MC-1 engine measurement variables 


n 


47 

164.5 

43 

505 

739.61 
765.54 
755.66 
790.43 
713.82 
167.86 
175.53 

141.29 

37.237 

875.68 

744.62 
725.24 
514.22 
65.915 

560.29 
81.133 

1607.7 

1442.7 

646.84 

l 18654 

t 48207 

m 


47 

164.5 

43 

505 

710.02 

736.12 

724.90 
758.69 
666.77 
167.57 
175.35 
136.49 

37.344 

833.91 
714.11 
696.15 
513.66 
63.227 

533.93 

77.972 

1533.4 

1385.4 

625.22 

i 18028 

46288 

ft 


47 

164.5 

43 

505 

750.52 

779.99 

767.42 

801.62 

722.84 

167.97 

175.51 

141.86 

37.102 

886.04 

759.20 

734.49 

514.28 

67.748 

568.59 

79.001 

1600.3 

1433.0 

649.86 

i 18839 

: 47301 

g§ 


47 

164.5 

43 

505 

725.22 

754.73 

744.82 

774.97 

683.53 
167.77 
175.44 
138.60 

37.303 

855.30 

736.26 

721.88 

513.86 

65.945 

549.53 
77.520 

1542.0 

1389.0 

627.10 

l 18390 

! 46142 

|§ 


47 

164.5 

43 

505 

725.22 

754.73 

744.82 

774.97 

683.53 
167.77 
175.44 
138.60 

37.303 

855.30 

736.26 

721.88 

513.86 

65.945 

549.53 
77.520 

1542.0 

1389.0 

627.10 

: 18390 

i 46142 

m 


47 

164.5 

43 

505 

735.96 

765.31 

756.07 

785.06 

700.38 
167.75 
175.35 
137.71 

36.953 

817.33 

761.46 

731.38 

513.46 
69.871 

554.00 
78.774 

1542.6 

1390.0 

645.62 

: 18382 

i 46696 

I 


47 

164.5 

43 

505 

735.96 

765.31 

756.07 

785.06 

700.38 
167.75 
175.35 
137.71 

36.953 

817.33 

761.46 

731.38 

513.46 
69.871 

554.00 
78.774 

1542.6 

1390.0 

645.62 

18382 

46696 

1 

Units 

psia 
deg R 
psia 
deg R 

psia 
psia 
psia 
psia 
psia 
deg R 
deg R 
Ib/s 

psia 
psia 
psia 
psia 
deg R 
Ib/s 

psia 
psia 
deg R 
deg R 

psia 

rpm 

lb 


Variable 

Name 

PSVL10 

TTVL10 

PRSPRV 

TTRPFV 

PSOXDS 

PSVL13 

PTVL14 

PSVL15 

PTVL18 

TTVL14 

TTVL18 

WOXTOTL 

PSVL00 

PSVL01 

PTVL05 

PTVL09 

TTVL05 

WRPTOTL 

PTHTGI 

PTVL22 

TTHTGI 

TTHTGD 

PTMCHY 

SNSHFT 

FT15A 

Subsystem 

Inlet 

LOX 

RP (Fuel) 

GG/Turbine 

MCC/Nozzle 

Other 


G-3 




Table G4 Baseline values for MC-1 engine hardware parameters 


m 


9.1375E-03 
9.8177E-04 
2.9281 E-05 
1.0155E-04 

1 .45768 
0.88478 
1 
1 

0.54068 

0.70188 

0.83975 

0.77519 

0.97501 
1 .04626 
1 
1 

0.78028 

0.03555 

1 

0.96998 
! 4.74871 

s 


9.1968E-03 

9.7986E-04 

3.1865E-05 

1.0610E-04 

1 .48744 

0.94156 

1 

1 

0.55772 

0.71219 

0.83940 

0.78249 

0.98886 

1.06951 

1 

1 

0.77630 

1 

0.03579 

1 

0.97559 
i 4.50382 

S 


9.0462E-03 
9.7464E-04 
2.8069E-05 
1 .0705E-04 

1 .50258 
0.88393 
1 
1 

0.55021 

0.70981 

0.81671 

0.74874 

0.97507 

1.03617 

1 

1 

0.81232 

1 

0.03901 

0.97675 

0.94017 

i 4.71635 

§§ 


9.0988E-03 
9.7336E-04 
2.8791 E-05 
1 .0242E-04 

1 .58894 
0.86265 
1 
1 

0.55176 

0.73173 

0.79169 

0.73115 

0.98478 
1 .04857 
1 
1 

0.80241 

1 

0.03108 

0.99648 
i 0.94960 

i 4.56995 

n 


9.0988E-03 
9.7336E-04 
2.8791 E-05 
1 .0242E-04 

1 .58894 
0.86265 
1 
1 

0.55176 

0.73173 

0.79169 

0.73115 

0.98478 
1 .04857 
1 
1 

0.80241 

1 

0.03108 

0.99648 

0.94960 

i 4.56995 

n 


9.2721 E-03 
9.7323E-04 
2.6701 E-05 
1.0123E-04 

1 .98884 
1 .04932 
1 
1 

0.54786 

0.70282 

0.81972 

0.74981 

0.96658 
1 .05700 
1 
1 

0.79622 

1 

0.03365 

0.97775 

0.94666 

i 4.57350 

1 


9.2721 E-03 
9.7323E-04 
2.6701 E-05 
1.0123E-04 

1 .98884 
1 .04932 
1 
1 

0.54786 

0.70282 

0.81972 

0.74981 

0.96658 
1 .05700 
1 
1 

0.79622 

1 

0.03365 

0.97775 

0.94666 

4.57350 

8 

Units 

lbf-s 2 /lbm-in 5 

lbf-s 2 /lbm-in 5 

lbf-s 2 /lbm-in 5 

lbf-s 2 /lbm-in 5 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

none 

Btu/sec 


Variable 

Name 

RCALMF 

RCALMO 

RKFL1 

ROLN1 

XMGGKO 

XMGGOO 

XMMCKO 

XMMCOO 

CDGGKI 

CDGGOI 

CDKINJ 

CDOINJ 

PSIMKPMP 

PSIMOPNIP 

TRQMKPMP 

TRQMOPMP 

CDGGNZ 

ETAMHTGT 

FRICFACT 

CDNOZL 

ECSMMCHB 

QDOTVL18 

Component 

Lines 

Valves 

Injectors 

Pumps 

GG/Turbine 

Turbopump 

MCC/Nozzle 

Other 


G-4 




Table G5 Revised hardware parameter list for MC-1 engine reduction study 



G-5 






G-6 



Appendix H 

MC-1 engine 

R2 and R3 test series data 


H 





H-l 



Table H2 MC-1 engine - test R2-2 one second average data 
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Table H3 MC-1 engine - test R2-3a one second average data 
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Table H4 MC-1 engine - test R2-3b one second average data 
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Table H5 MC-1 engine - test R2-4 one second average data (continued) 
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Table H5 MC-1 engine - test R2-4 one second average data (continued) 
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Table H5 MC-1 engine - test R2-4 one second average data (continued) 
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Table H5 MC-1 engine - test R2-4 one second average data (continued) 
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Table H5 MC-1 engine - test R2-4 one second average data (continued) 
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Table H6 MC-1 engine - test R3-1a one second average data 
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Table H8 MC-1 Engine - R2 series tests - general observations 
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Table H9 MC-1 Engine - R2 series tests - observations adjusted to std input 
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Appendix I 

MC-1 engine 

GDR prediction comparisons 
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Table 12 Comparison of GDRA (linear) and GDRC(2nd order) results 



1-2 


% Dev - average percent absolute deviation from standard ROCETS data reduction results 



Table 13 Comparison of GDRA and ROCETS/DR results using modified hardware set 
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% Dev - average percent absolute deviation from standard ROCETS data reduction results 




Table 14 Sensor elimination study results 
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% Dev - average percent absolute deviation from standard ROCETS data reduction results 



























Table 14 Sensor elimination study results (continued) 
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Table 15 Comparison of GDRA and ROCETS/DR results using flight measurements only 
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Table 16 Subset selection examples with MC-1 engine parameters 
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